Spin Correlations in Ho 2 Ti 2 07: a Dipolar Spin Ice System 
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The pyrochlore material Ho2Ti2 07 has been suggested to show "spin ice" behaviour. We present 
neutron scattering and specific heat results that establish unambiguously that Ho2Ti2C>7 exhibits 
spin ice correlations at low temperature. Diffuse magnetic neutron scattering from Ho2Ti207 is 
found to be quite well described by a nearest neighbour spin ice model and very accurately described 
by a dipolar spin ice model. The heat capacity is well accounted for by the sum of a dipolar spin ice 
contribution and an expected nuclear spin contribution, known to exist in other Ho 3+ salts. These 
results settle the question of the nature of the low temperature spin correlations in Ho2Ti2 07 for 
which contradictory claims have been made. 



Spin ice materials |l],|^] are magnetic substances in 
which the atomic magnetic moments obey the same or- 
dering rules as the hydrogen atoms in ice, H2O. They 
provide a bridge between the simple statistical mechan- 
ics of ice- type models [ | , and the complex behavior of 
frustrated magnets with wider relevance to diverse 

areas of research such as high temperature superconduc- 
tivity and neural networks. Furthermore, relatives of spin 
ice such as R2M02O7 (R = rare earth) are a current fo- 
cus of attention for their interesting electronic properties 
flpl. It i s therefore important and desirable to estab- 
lish the detailed physics of the simplest spin ice materi- 
als, which provide among the best model systems for the 
study of frustration and its broad consequences. 

Experiments on Ho2Ti207 provided the original moti- 
vation for the concept of spin ice 0]. In cubic Bx>2Ti 2 07, 
(space group FdZm), the magnetic Ho 3+ ions occupy a 
cubic pyrochlore lattice, a corner-linked array of tetra- 
hedra that is identical to the lattice formed by the mid- 
points of the oxygen-oxygen bonds of cubic ice Jl(| . The 
ground state of Ho 3+ is an effective Ising doublet of lo- 
cal (111) quantization axis []T],p~l] p~4| , with nearest neigh- 
bor Ho 3+ moments experiencing an overall ferromagnetic 
coupling that is largely dipolar in origin. Qualitatively, 
Ho2Ti207 thus approximates the nearest neighbor spin 
ice model ]l|, p^Ul6| , in which ferromagnetically-coupled 
Ising spins on the pyrochlore lattice lie parallel to the lo- 
cal (111) axes that point towards the centers of the tetra- 
hedra. The ground state of this model requires two spins 
pointing into and two pointing out of each tetrahedron. 
Then if spins represent hydrogen atom displacement vec- 
tors, one obtains the ice rules that lead to Pauling's result 
for the extensive zero temperature entropy of the disor- 
dered ground state of ice 0Jl^| • 

Neutron scattering experiments on Ho2Ti207 in zero 



magnetic field show only broad diffuse scattering down to 
a temperature T ~ 0.35 K, consistent with a disordered 
spin ice state [|jl8 , while low temperature muon spin 
relaxation (/iSR) work also finds no evidence for a mag- 
netic transition Field- induced ordered states are 
also consistent with the spin ice scenario |lj|ll|. Com- 
bined with a ferromagnetic Curie- Weiss temperature of 
dew — 1.9 ± 0.1 K, Ho 2 Ti 2 07 would thus appear to 
be a prototypical spin ice material. Another candidate 
for spin ice behaviour is Dy2Ti207 ||, for which spe- 
cific heat measurements down to 250 mK can be inte- 
grated to give the residual entropy expected for the spin 
ice model. However, none of the measurements reported 
in Ref. |lj or Ref. Q constitute a proof that the spin 
ice ground state exists in these materials in zero ap- 
plied field. Moreover, a certain amount of confusion has 
recently arisen due to some contradictory observations: 
Siddharthan et al. suggested that the specific heat be- 
haviour of Ho2Ti207 possibly indicates a transition to a 
partially ordered state at ~ 0.8 K [jl9|, in disagreement 
with the neutron scattering and /iSR results Jl|Jl^|. In 
this Letter we present neutron scattering and new heat 
capacity data that unambiguously establish the spin ice 
nature of the zero field spin correlations in Ho2Ti207. 

The diffuse magnetic neutron scattering from a flux 
grown single crystal of Ho2Ti207 was measured in the 
static approximation on the PRISMA spectrometer at 
ISIS. The crystal was oriented with [1,-1,0] vertical such 
that the (h, h, I) scattering plane included the three prin- 
cipal symmetry axes [1,0,0], [1,1,0] and [1,1,1]. Figure 
la shows the scattering pattern at T ~ 50 mK. One of 
the main features of the experimental data is the 'four- 
leaf clover' of intense scattering around (0,0,0). There 
is also strong scattering around (0,0,3) and a broad re- 
gion of slightly weaker scattering around (3/2,3/2,3/2). 
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These intense regions are connected by narrow necks of 
intensity giving the appearance of bow-ties. The width 
of the intense regions indicates short range correlations 
on the order of one lattice spacing. Qualitatively similar 
scattering has been observed in ice itself Q . 

To model this data we use the standard expression for 
the neutron scattered intensity I(q) |2l|], along with the 
Hamiltonian It 
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where Ising spins of unit length are constrained to 
their local z\ = (1,1,1) axes; J is a near neighbor ex- 
change coupling and D the dipolar coupling. Because of 
the local Ising axes the effective nearest neighbor energy 
scales are J nn = J/3 and D nn = 5D/3 |24|. 

The near neighbor spin ice model |l|,|l6| corresponds 
to D = and J positive (ferromagnetic). Data were 
simulated with 3456 spins (6x6x6 cubic unit cells) at 
T/J = 0.15 where the model is effectively in an ice-rules 
ground state. The calculated pattern is shown in Fig. lb. 
It successfully reproduces the main features of the exper- 
imental pattern, but there are differences, notably in the 
extension of the (0, 0, 0) intense region along (h, h, h) and 
the relative intensities of the regions around (0, 0, 3) and 
(3/2,3/2,3/2). Also, the experimental data shows much 
broader regions of scattering along the diagonal direc- 
tions. Clearly, the experimental spin correlations do not 
reflect a completely disordered arrangement of ice-rule 
states, but some states are favored over others p^,p^|. 

The more complete dipolar spin ice model p2|-|24[, 
has D nn = 2.35 K, fixed by the lattice constant, and 
Jnn = —0.52 K a negative (antiferromagnetic) parameter 
determined by fitting the peak temperature of the elec- 
tronic magnetic heat capacity (see below). Spin ice be- 
havior emerges in this model from the dominant effect of 
the long-range nature of dipolar interactions [p2T-p4f , and 
can quantitatively account for the heat capacity data of 
Dy2Ti2C>7 (taken from Ref. 0) |22]fp|l. Using single spin 
flip dynamics, J(q) was calculated |21| on a system size 
of 1024 spins (4 x 4 x 4), at T = 0.6 K where significant 
ground state correlations have developed, using standard 
Ewald summation techniques which properly handle in- 
finite summation of dipole-dipole energy terms (see Refs. 
p^HJ). The calculated pattern is shown in Fig lc. It 
captures most details of the experimental pattern missed 
by the near neighbor spin ice model in Fig. lb, such as 
the four intense regions around (0,0,0), the relative in- 
tensities of the regions around (0, 0, 3) and (3/2, 3/2, 3/2) 
and the spread of the broad features along the diago- 
nal. Note also the low scattering intensity at (2,2,0), 
consistent with the experimental low intensity around at 
(2, 2, 0). The approach of J(q) to zero at q = is consis- 
tent with a recent observation that the ac-susceptibility 
of Ho2Ti2C>7 approaches zero below 0.7 K p5||. 
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FIG. 1. (a) Experimental neutron scattering pattern of 
Ho2Ti207 in the (hhl) plane of reciprocal space at T ~ 50 
mK. Dark blue shows the lowest intensity level, red-brown the 
highest. Temperature dependent measurements have shown 
that the sharp diffraction spots in the experimental pattern 
are nuclear Bragg peaks with no magnetic component, (b) 
J(q) for the nearest neighbor spin ice model at T — 0.15J. 
(c) /(q) for the dipolar spin ice model at T — 0.6 K. The 
areas defined by the solid lines denote the experimental data 
region of (a). 
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The difference between Figs, fb and fc shows that 
dipolar interactions, whilst inducing a low energy ice- 
rules manifold, do cause further correlations to emerge 
among the spins than those that are due to a nearest 
neighbor ferromagnetic interaction. This thermal bias to- 
wards certain ice-rule configurations in the dipolar model 
is consistent with mean field calculations ^4|, and also 
with the behavior observed when single spin flip ther- 
mal barriers are removed by using multiple-spin flip al- 
gorithms in numerical simulations |25j . In that case, 
the intense regions of I(q) evolve into Bragg peaks upon 
cooling below 0.17 K as the system enters a long range 
ordered state. It may be that at low temperature in 
Ho2Ti2C>7, short range correlations associated with a pu- 
tative ordered state begin to emerge, but are dynamically 
inhibited from developing into Bragg peaks upon further 
cooling. 

To further test the dipolar model, we show in Fig. 2 
a comparison between the calculated I(q) and that mea- 
sured in a separate experiment with the same crystal on 
the IN14 spectrometer at the ILL, Grenoble. Here it was 
found that the diffuse scattering was nearly temperature- 
independent below 0.7 K. For the theoretical fit, the form 
factor was determined empirically by fitting the neutron 
data at 35 K in the paramagnetic regime . The com- 
parison between theory and experiment obtained at 0.7 
K along (0, 0, 1) is shown in Fig. 2a, where it is seen 
to be very satisfactory. Hence it is clear that the dipo- 
lar Hamiltonian (1) describes the magnetic behaviour of 
Ho 2 Ti 2 07 to a very good approximation. 
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FIG. 2. Experimental neutron scattering intensity, /(q), of 
Ho2Ti2C>7 (filled symbols) along the (001) direction of recipro- 
cal space. For quantitative comparison is the intensity (open 
symbols) obtained from Monte Carlo simulation of the dipolar 
spin ice model with J nn = —0.52 K and D nn = 2.35 K. 

We now turn to the heat capacity, analysis of which is 
complicated by the previously reported difficulty in equi- 
librating this material at low temperatures (T < 0.8 K) 
p9[ . Single crystals of Ho2Ti2C>7 were synthesized using 
the floating-zone method described in Heat capac- 
ity measurements were performed in a Quantum Design 
PPMS system inside a 9T superconducting magnet. The 
heat capacity was measured using a thermal relaxation 



method between 0.34 and 20 K. Good thermal contact 
between the single crystal and sample stage was ensured 
with the uses of low temperature grease. Our heat capac- 
ity data, shown in Fig. ^, go to lower temperature than 
those of Ref . jl9) and show no evidence for a phase transi- 
tion. Rather, there is a noticeable change in dynamics as 
equilibration times become longer deep within the spin 
ice regime, as signaled by the vanishing electronic spe- 
cific heat capacity, but the total heat capacity continues 
to steadily increase as the temperature is lowered (note, 
however, the shoulder at T ~ 1.5K). 

Blote et al. measured the heat capacity of the py- 
rochlore Ho 2 GaSb07 [llj and obtained similar behav- 
ior to that found for Ho2Ti 2 07. They successfully ac- 
counted for the heat capacity below ~ 2 K by introduc- 
ing a Schottky anomaly with a theoretical maximum for 
Ho (/ = 7/2) of 0.9 R, due to a splitting of the eight 
nuclear levels with a level spacing of 0.3 K. 
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FIG. 3. The total specific heat of Ho2Ti207 is shown by the 
empty squares. Expected nuclear contribution is indicated by 
the line, while the electronic contribution has been estimated 
by subtracting these two curves (open circles). Near 0.7 K 
this subtraction is prone to a large error (see text). Dipolar 
spin ice simulation results are indicated by the filled circles. 

The electronic component of the magnetic heat capac- 
ity for Ho2Ti207 may be extracted by subtracting off the 
nuclear Schottky contribution estimated by Blote et al. 
for Ho2GaSb07 |11|. This reveals an electronic magnetic 



heat capacity that shows the same characteristic spin ice 
shape as observed for Dy2Ti207 [||. Fig. || shows a very 
broad electronic magnetic peak centered at T ~ 1.9 K. 
With this heat capacity revealed, a natural explanation 
of the observed change in dynamics in this material at 
T ~ 0.8 K, (attributed to a transition in Refs. (li|2f]), is 
that this is approximately the temperature below which 
macroscopic relaxation of the magnetic degrees of free- 
dom into the low energy ice-rule manifold has essentially 
ceased (this is also consistent with Rcf.( [p5[)). In other 
words, the onset of spin ice physics is much higher than 
this temperature, most likely T ~ 1.9 K. The behavior of 
Ho2Ti207 may be compared to that of Dy2Ti207, where 
the onset of spin ice correlations (peak in the specific 
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heat) is ~ 1.2 K and the magnetic heat capacity is al- 
most zero at T ~ 0.25 K g. 

The subtraction of the nuclear contribution, itself an 
approximation, becomes prone to error in the range 0.6- 
0.9 K where the nuclear specific heat rises sharply. Hence, 
although our subtraction reveals a small feature in the 
electronic specific heat at T ~ 0.7 K, a subtraction of the 
nuclear contribution from the total heat capacity data 
reported in Ref. gives no such feature; thus we con- 
clude that an electronic-only heat capacity in this region 
cannot be reliably determined. 

Based on the electronic specific heat peak height, and 
the temperature at which it occurs, the dipolar spin ice 
model allows two independent methods for determining 
the value of the nearest neighbor exchange for Ho2Ti2 0y 
p2|. Using these procedures we find a nearest neighbor 
exchange J nn ~ — 0.52 K. In Fig. || we show specific 
heat results from a standard Monte Carlo single spin flip 
simulation of the dipolar spin ice model, for a system 
size of 6 x 6 x 6 unit cells (3456 spins) with interaction 
parameters J nn = —0.52 K and nearest neighbor dipole 
strength D nn = 2.35 K. We find quantitative agreement 
between the experimental electronic specific heat, C, and 
that found from the simulation. At high temperature, C 
approaches the dipolar spin ice heat capacity smoothly, 
as it should since the nuclear contribution is negligible. 
At low T it follows the nuclear contribution once the 
electronic component has died. Finally, in-between the 
low T and high T regimes, there is the shoulder at T ~ 
1.5 K in the total C(T) (alluded to above) corresponding 
to a shift of regime between the dipolar spin ice controlled 
regime and the nuclear hyperfine one. 

Taken with our neutron scattering results, it appears 
clear that Ho2Ti207 is quantitatively well characterized 
by the dipolar spin ice model with a missing entropy close 
to Pauling's prediction, and what is found for Dy2Ti2C>7 
P,p2j. The spin ice behavior revealed explicitly here 
for Ho2Ti207 is in contrast to the conclusions of Refs. 

These were based on a finite size truncation 
of the dipolar interaction in the simulations and only 
a qualitative comparison with experimental specific heat 
that neglected the important nuclear contribution (Fig. 
3 of ]19[| ). Our deduced value of the antiferromagnetic 
exchange interaction indicates that Ho2Ti207 is further 
into the dipolar spin ice region of the magnetic phase di- 
agra m p2f than is Dy2Ti2C>7, contrary to the claims of 
Ref. [[28|, with J^°/D^° ~ -0.22 and J^i/D^l ~ -0.53. 

In conclusion, we have unambiguously established that 
Ho2Ti207 possesses a spin ice state in zero field. Experi- 
ments investigating the spin correlations of Dy2Ti207 are 
in progress p9| . The Ho salt, in contrast to its Dy ana- 
logue, contains a single rare earth isotope that lends itself 
well to neutron scattering and the study of hyperfine ef- 
fects. With intriguing spin dynamics and field-induced 
ordering phenomena, Ho2Ti207 offers much new physics 
to be explored in the field of frustrated magnetism. 
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